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Abstract 

Because the source term for the equations of motion of a conformahy coupled scalar field, such as 
the dilaton, is given by the trace of the matter energy momentum tensor, it is commonly assumed 
to vanish during the radiation dominated epoch in the early universe. As a consequence, such 
fields are generally frozen in the early universe. Here we compute the finite temperature radiative 
correction to the source term and discuss its consequences on the evolution of such fields in the 
early universe. We discuss in particular, the case of scalar tensor theories of gravity which have 
general relativity as an attractor solution. We show that, in some cases, the universe can experience 
an early phase of contraction, followed by a non-singular bounce, and standard expansion. This 
can have interesting consequences for the abundance of thermal relics; for instance, it can provide 
a solution to the gravitino problem. We conclude by discussing the possible consequences of the 
quantum corrections to the evolution of the dilaton. 

PACS numbers: 
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I. INTRODUCTION 



Although fundamental scalar fields have yet to be discovered, their role in fundamental 
physics is unquestionable. One particular class of scalar fields are those which are confor- 
mally coupled to matter. Well known examples of such fields are the Jordan-Fierz-Brans- 
Dicke (JFBD) scalar which is introduced as an extension to general relativity (GR) [l| or 
the dilaton which in the context of string theory arises as part of the gravitational multiplet 
[2I. Of course string theories possess many other scalars or moduli which are conformally 
coupled to matter as well. 

Through a series of field redefinitions, theories with conformally coupled scalar fields 
(CCSFs) can always be reexpressed in terms of Einstein gravity with well determined cou- 
plings to the matter sector of the theory. These new scalar interactions may be perceived 
as new forces which can lead to violations of the equivalence principle if all couplings are 
not universal j^. Placed in a cosmological context, the evolution of these fields can lead 
to the variation of fundamental constants, including gauge couplings and Yukawa couplings 
in the general case. Theories with universal couplings to matter [3, 0] as well as so-called 
chameleon theories both of which are based on CCSFs, are constructed to minimize 
these violations (in particular they ensure the universality of free fall and the constancy of 
all non-gravitational constants). Indeed, precision tests of gravity often lead to important 
constraints on the parameters of the theory with these scalar fields. Cosmologically, the 
presence of a new dynamical degrees of freedom can also have severe consequences on big 
bang nucleosynthesis (BEN) 0, S 0, [13, 11, 12|. Thus there are a whole host of observa- 



tions, including cosmic microwave background anisotropies [IJ] and weak-lensing [15| , which 
constrain theories with CCSFs. 

The cosmological evolution of CCSFs are determined by the standard equation of motion 
for a scalar field with an extra source term proportional to the trace of the energy momentum 
tensor. That is, in the Einstein frame, we can write (see Section IIIII below for details and 
definitions) 

+ 3H^, = -^^^ - ^7iGM^*K (1) 
where star denotes quantities in the Einstein frame. 

, , , cilnA , , 

characterizes the strength of the scalar interaction and A{(j)^) is the coupling function, that 
appears in the conformal transformation relating the string or Jordan frame to the Einstein 
frame 

gf„y = A^{(j)^)g^^^. (3) 

One can recover Einstein relativity in a simple way by introducing a potential V that 
generates a sufficiently high mass for 0=,, and forces 0* to a constant value (note, however, 
that this value will depend on the local energy density as in the chameleon mechanism [^). 
Even in the absence of a potential, this theory will have GR as an attractor if there is a well 
defined minimum to In A In a JFBD theory, A oc e"'^"'^* and there is no attraction to GR, 
as in the case of the dilaton in heterotic string theory. In these models, a is constant and 
the deviation from GR is fixed (in terms of the PPN parameters). In this case, the theory 
is only compatible with GR if A is small enough. It is also clear from Eq.([T]) that during the 
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radiation dominated epoch in the early universe, the field will remain constant until 7^ 0, 
which normally would correspond to the onset of matter domination 1^. This effect has 



been utilized, for example, as a potential solution to the problem of a run away dilaton jjj 

Assuming a perfect fluid form for the energy momentum tensor, its trace is simply, p — 3p, 
and vanishes in a radiation dominated epoch when the equation of state is characterized by 
p = p/3. While this is certainly a good approximation deep in the radiation era, it breaks 
down whenever, the ternperatur e ap proaches a mass threshold, and some particle species 
become non-relativistic j3, 0, [O, 12]. There, the trace becomes proportional to 



^2 ^2 



where g is the number of degrees of freedom for the particle at threshold, the sign refers 
to different statistics (+ for fermions, — for bosons), T is the temperature of the radiation 
in the Jordan frame, and z = m/T. The source terms from particle thresholds and in 
particular that of the electron, was shown to greatly relax the bounds from BBN on the 



initial conditions of theories of gravity with CCSFs ^2 . 

A non- vanishing source term also arises when one includes the contribution from the trace 
anomaly [6,[l3|. Although this is generally small - we show that w — 1/3 = O (10~^ — 10~^) - 
there can be situations in which this contribution has a signiflcant impact on the cosmology 
of these models. In fact, not only can it affect the evolution of the scalar fleld, but we 
show that, in some cases, it can lead to a brief early phase of contraction followed by a 
non-singular bounce in the cosmological scale factor. In these cases, one can also derive a 
maximum temperature of the universe which may provide a solution to the gravitino problem 
or other unwanted relics. 

The paper will be organized as follows. In the next Section, we will flrst calculate the 
contribution to the equation of motion from the trace anomaly. This will be done in both 
the standard model (SM) and the minimal supersymmetric standard model (MSSM). We 
will then set up our formalism for treating gravity theories with CCSFs in Section 3. In 
Section 4, we will present the results for the cosmological evolution of the scale factor and 
scalar fleld. A discussion and conclusions will be given in Section 5. 



II. EQUATION OF STATE OF THE SM AND THE MSSM AT HIGH TEMPER- 
ATURES 

From Eq. ([1]), the source term for the equation of motion for a CCSF, is proportional 
to p — 3p = p{l — 3w) where the parameter w characterizes the equation of state. In the 
radiation dominated era, and in the absence of interactions, w = 1/3. Here, we compute 
corrections to w in the context of the SM as well as the MSSM at temperatures higher than 
the masses of the particles. This computation is readily extended to scalar-tensor theories of 
gravity, since they satisfy the weak equivalence principle (namely, non-gravitational physics 
is unaffected by the scalar fleld). This is manifest in the Jordan frame, in which matter is 
not directly coupled to the scalar. However, it is also easy to interpret the computation in 
the Einstein frame [see the discussion after Eq. ([SI])]. It is useful to start from the expression 
for the free energy density, which we parametrize as 

^=-3/gT'. (5) 
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In this expression, gf is the effective free energy number of relativistic degrees of freedom, 
that is defined by normahzing to the free energy density of one non-interacting massless 
bosonic degrees of freedom: 

^l,free = (6) 

For vanishing chemical potential, the pressure and energy density can be obtained from 
Eq. © as 



From these expressions, one finds 



P-3P = -T^#^, (8) 
^ 90 dlnT ' ^ ^ 

or, equivalent ly, 

P ( dlngfY^ 

We see that the trace of the energy-momentum tensor ([8]) is different from zero when- 
ever the effective number of degrees of freedom changes with the temperature. Most early 
universe calculations in cosmology consider only the contribution of the relativistic species 
to the pressure and energy density (since those of massive species are exponentially sup- 
pressed). Therefore, one includes only particles with mass smaller than the temperature in 
the expressions ([7]), and, as a consequence, gf only varies whenever the temperature drops 
below any mass threshold. However, g^ also varies at temperatures above all mass thresholds 
in the theory, due to the trace anomaly j^. To compute this, we recall the relation 

T=~\nZ (10) 

where Z is the partition function and V the physical volume (one typically imposes periodic 
boundary conditions on a cube of total volume V\ assuming homogeneity on a sufficiently 
large volume, the pressure and energy density do not depend on V). The partition function 
can then be computed diagrammatically [181]. The one loop vacuum diagram provides the 
free theory result (jH]). Once we factor out the numerical factor — 7r^/90, the coefficient gf is 
simply the number of bosonic and fermionic relativistic degrees of freedom 

^;,free = N, + N f , (11) 

where the 7/8 coefficient multiplying the fermionic contribution is due to the different value 
of the partition function obtained from the Fermi-Dirac thermal distribution rather than 
the Bose-Einstein distribution. In the SM, we have N\, = 28 bosons and Nf = 90 fermions, 
giving ^/,free = 427/4 = 106.75. 

Higher loops account for the interactions. The departure from the ideal gas regime due to 
the interactions reduces the effective number of degrees of freedom. Here, we only consider 
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the leading order correction encoded in two loop vacuum diagrams.^ The corrections are 
proportional to the gauge, Yukawa, and Higgs couplings. The variations of these couplings 
with temperature then lead to the nontrivial dependence of (?/ on the temperature, and, as 
a consequence, to p — 3 P 7^ (see Fig. [T]). 




a3(MSSM) 




FIG. 1: Left panel: The absolute value of the leading corrections (two-loops) to the effective 

free energy number of relativistic degrees of freedom from the different gauge interactions, top 

Yukawa and Higgs self-coupling in the SM. The strong contribution (03), is dominant and is also 

plotted for the case of the MSSM assuming that all supersymmetric partners have masses fixed 

(2) 

at msusY = 500 GeV. All contributions to are negative. Right panel: The absolute value of 

the leading correction to 1/3 — uj from the same interactions shown in the left panel. The strong 

(2) 

contribution is dominant not only due to the fact that the correction to is bigger, but also 
because the running of 03 is steeper. On the other hand, the contribution from the hypercharge 
interaction (and the Higgs coupling for high temperatures) is negative and opposite to the others 
since ai (A) increases with T. 



For the SM, the one and two loop contributions give [18|, ll9|, l20|, l2l|, l22|, l23|, [24] : 

9f 
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(12) 



87r2 



where 



32 7r2 



I7/ P 



VapV (3a 1 



^ As we shall see, the effect is dominated by the strong interactions; at temperatures just above the masses 
of the SM particles, higher order loops are also important. The various loop contributions for QCD, 
up to aa"^ ln(l/Q;3) have been evaluated in Ref. 19 1, where it is shown that only for temperatures T > 
lO^GeV, the two loop contribution dominates higher order terms. As we are interested in the dynamics 
of the Universe at very high energies, we will discuss only the leading term and assume that is a good 
approximation. 
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are the Yukawa matrices of the SM; A is the quartic Higgs couphng, and as, 02, cty are the 
couphng constants associated with the strong, weak and hypercharge interactions respec- 
tively. Our results are expressed in terms of the coupling ai = 5ay/3, (we are following 



notation from Ref. [25|). Explicitly: 



%y 5a a gl . 

IzTT 3 cos^ Ow 4:71 sm 9w 47r 

where a = e^/47r, and gy, g2 and 5^3 are the usual Uy{1), SU{2) and SU{3) coupling 
constants. We are using the modified minimal subtraction {MS) scheme [26i] as a renormal- 
ization prescription. 

Results from LEP and the Tevatron imply [27|: 



ai{Mz) - 0.017, a2iMz) ~ 0.034, asiMz) ~ 0.118. (14) 

At unification scales (T ~ lO^^GeV), the electroweak couplings are comparable to the strong 
coupling and their contributions could be more important. Eq. f[T^ and Fig. [T] show, 
however, that this is not the case, and the strong contribution (03) is dominant even at very 
high energies due to the large number of strongly interacting particles. In addition, the third 
generation dominates the contribution of the Yukawa couplings. By taking into account the 



SM masses |27|, 1281]: mt ~ 173 GeV, nib ^ 4.2 GeV and ~ 1.78 GeV, it is evident that 



we can neglect all the Yukawa contributions apart from the top: 

2 

|l/^^(m,)|2^2^^ 0.99, (15) 

V 

where the vacuum expectation value of the Higgs is w ~ 246 GeV. Finally, A is the only 
unknown parameter of the SM. In the simplest model. It is related directly to the Higgs mass: 
A = (m/i/f)^. Combined results from the four LEP collaborations (ALEPH, DELPHI, L3 
and OPAL) gives the lower bound [23]: mh > 114.4 GeV at 95% of confident level, which 
implies 

A(m,,) > 0.216. (16) 

In Fig. dl we have fixed rrih = 120 GeV. 

As mentioned earlier, the temperature dependence of the couplings leads to the deviation 
from an ideal relativistic gas. This dependence is described by the renormalization group 
equations (RGE). Taking the leading (one- loop) contribution |25| to the running, we find: 

= -J— (-S423a^^ -135 aia2 + 9875 al + 58800 afi 
amT 80 vr^ ^ ^ 



15 /_ _ . _ 3 2 \ 2 



17ai + 45 a2 + 160 03 + — yt y 
256 TT'^ V ZTv / 

9A / 5 \ 

The contribution to dgj/dlnT from the strong interaction is even more dominant than the 
contribution to gf, since the running of is also more important than that of the other 
couplings. In this case, keeping only the QCD interaction is a good approximation up 
to order 10%, over the entire range of energies considered, unless the Higgs is very heavy 
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{nih ^ 170 GeV). At high energies the quartic Higgs interaction (and the hypercharge 
couphng) runs in the opposite direction to that of as, which partially compensates for the 
QCD contribution. 

The greatest source of uncertainty in this computation is the effect of higher order QCD 
terms in Eq. (IT^ . To estimate this effect, we included terms up to O (as*^^) in the com- 
putation of Qf for the SM (the various terms are listed in Ref. [ll]), and computed the 
corresponding value for w. We found a value for 1/3 — w about 50% smaller than that 
obtained from the O (0:3) term alone reported in f|T2l) . in the entire range of temperatures 
10 GeV < T < 10^^ GeV. Unfortunately, the next order contribution has an undetermined 
coefficient, and it has been shown in Ref. [l^ that this could have a non negligible effect on 
Qf. Moreover, there are nonperturbative contributions which could also be substantial. Due 
to these uncertainties, we choose to retain only the first nontrivial, O (as), contribution to 

"comforted" by the fact that, summing up to the highest fully known order, gives a result 
which is in reasonable agreement with the one we adopt. 

In summary, for rrih ^ 1 70 GeV, we have the approximate result 

= - 420 as + 0(af) , ^ = 2940 a^ + O (af) , 

1 _ 560 ^2 , ~ rrr\ - '^3 ^, ^3 {m) 



3 183 V ^ 27r i + 7tt3(m0 ln(^J 

with as (mt) ~ 0.0172. 

We can repeat this computation for the MSSM. Including only the (super)-QCD contri- 
bution, we find 

^ _ 1 son ^„ ^ n r^3/2\ dgj . , r^5/2\ 



^/ = — - 1890as + O(af ) , ^ = 5670a^ + O(a 
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w — 



1 168 as , ^o/-5/2^ ~ .^.x «3 ^SUSY .^Qv 

Q= ^ + ^l°^3 J , as(T)~- — ^ / ^ N , 19 

3 61 V / l + Sas msusY ln( ^ ^ 



msuSY 



Assuming for definiteness that the supersymmetric partners have masses msusY = 500 GeV, 
and using the running in Eq. ( fTSl) from T = niz toT = msusY, we have as(msusY) — 0.0153. 

We recall that the expressions f|T8|) and f|T9l) hold for temperatures greater than the 
masses of the particles, i.e. for T 200 GeV in the SM, and (with our assumption) for 
T >j 500 GeV in the MSSM. In the right panel of Figure [1] we show different contributions 
to the departure of the equation of state from the value w = 1/3 value corresponding to 
the noninteracting case. We note that the dominant effect decreases as the temperature 
increases, since as decreases with temperature (and, consequently the absolute value of the 
two loop contribution to gf decreases). We also notice that, at the lowest temperatures 
shown, the departure is stronger in the SM than in the MSSM case. This is because as runs 
faster in the SM . However, since as decreases less in the MSSM, the departure from the 
noninteracting value becomes greater for the MSSM as the temperature increases. 



III. SCALAR-TENSOR THEORIES OF GRAVITY 

In scalar-tensor theories, gravity is mediated by both a spin-2 graviton and a spin- scalar 
field that couples universally to matter. Following the notation of Refs. H, 14, 15 1, we start 
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from the action 
S 



-9 



R 



g^''d^<Pd,(P-2U (0) 



(20) 



where the first term is the action for the spin-2 graviton and the scalar field, while the second 
term is the action for matter. This expression defines the theory in the Jordan frame, in 
which the standard expression of the metric g^i, is used in the action for the matter fields 
(where by "matter" we generally denote any field in the theory apart from and gfj,u)- 
However, the gravitational interaction between matter fields is modified, since the function 
A {(f)) multiplies the Ricci scalar. Such a theory depends on 2 arbitrary functions and we 
shall assume in this work a vanishing potential for the scalar field, U = . We emphasize 
that is the bare gravitational constant and does not correspond to the gravitational 
constant that would be measured in a Cavendish- type of experiment^. 
Under a conformal transformation 



= A ^ (0) g^^ , 

and for a vanishing potential, the action of the system becomes 

d^x 



S 



where the scalar field 0* is defined by 



A' (0.) g, 



fii/ 



(21) 



(22) 



'd\nA{(f))Y ^2(0) 



(23) 



and A (0^) is short for A (0 (0*)) . 

The expression fl22l) is the action of the system in the Einstein frame. It is characterized 
by a standard action for the spin-2 graviton; however, the combination A^ (0*) g^fj,^, rather 
than the metric itself, is used in the action for matter. While the two expressions are 
equivalent, the Einstein frame is more often used to study the cosmology of the system 
(since the resulting cosmological equations are the standard ones), while the Jordan frame 
is more often used to study particle physics processes (since the physical lengths and masses 
are constant in this frame). Moreover, the Einstein frame is built in such a way that the 
kinetic term of the spin-2 and spin-0 mediators is diagonal so that the Cauchy problem is 
well-paused in this frame. 

We consider a FLRW universe with Euclidean spatial sections, in the Einstein frame, the 
line element is 

ds'^ = -dtl + Rl (t) dx' dx' (24) 



^ It can be shown that Newton's constant is given by Gn — G*^^(0*)(1 + a^) which depends a priori on 
time through 0* . The deviations from general relativity can be evaluated in terms of the post- Newtonian 
parameters 



as long as the field is light enough, as considered in this work; (see Ref. [2 
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and we denote by p* and P*, respectively, the total energy density and pressure of the matter 
fields. We then find the evolution equations |12| . 



3H^ = SttG.p, + 



3 (PRif 



47rG, (p, + 3P,) + 2 



R* dtl \ dt^ 



+ 3H^-— = -A-rvG^a (0*) (p* - 3P* 



dt^^ *dt, 

^ + 3H, (p, + P.) = a (0.) (p. - 3P,) ^ (25) 

where the Hubble rate is defined as H^, = rflnP^/rft*, and a (0*) was defined in Eq. ([2]). The 
first Eq. of fl25l) is the 00 Einstein equation while the second equation is a linear combination 
of the 00 and ii Einstein equations. We note that these equations are the standard FLRW 
equations, and the scalar field contributes to them only through its energy density and 
pressure. The other two equations are the field equations for the scalar field, and for matter, 
respectively. One out of the last three equations can be derived from the others in Eq. (I25p . 
as a consequence of a nontrivial Bianchi identity. Both source terms are proportional to a 
and refiect the scalar interaction and the action-reaction law (or the fact that the sum of the 
two stress-energy tensors is covariantly conserved). In the third equation, we see that the 
scalar field is coupled to the trace of the energy momentum tensor of matter as described 
in the introduction. This term is usually assumed to vanish in the radiation dominated era. 
However, as we discussed in the previous Section, this term is actually always nonvanishing 
at some level. Finally, we see from the last expression in Eq. (!25l) that the scalar field also 
modifies the continuity equation for matter. 

The geometry in the Jordan frame is also of the FLRW type, with scale factor and time 
related to those in the Einstein frame by 

R = A{(t),)R, , dt = A dt,. (26) 

As a "time variable" for the evolution, we choose the e-folds of expansion p in the Einstein 
frame (not to be confused with the pressure P). Denoting by P*,in the value of the scale 
factor at some initial time (for instance, at the end of reheating, when the thermal bath has 
formed), the number of e-folds is 

P - . (27) 

From the relations ( l26l) and ( 1271) . one can immediately relate the derivatives with respect to 
the time in either frame to derivatives with respect to p: 

d 1 d A d , , 
(IR^ 

dp dt^: dt 

In what follows, prime will denote differentiation with respect to p . 

It is worth noting that P* grows monotonically with time, since is governed by the 
standard Friedmann equation (the first of Eq. (B3i) ) and it is strictly positive, hence ensuring 

does not change sign so that p{t^) is a monotonic function. Therefore, we can indeed use 
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p as "time variable" . On the contrary, the Hubble rate H in the Jordan frame is related to 
that of the Einstein frame by 

H=^{l + a<j>:). (29) 

As we shall see in the next Section, it is possible that H becomes negative, indicating 
that the universe is actually contracting in the Jordan frame for some time during the 
evolution. Therefore, the scale factor and the temperature in the Jordan frame can have a 
non monotonic evolution, and we cannot use either as "time variables" for the system. 

The evolution equations in the Jordan frame are given in Ref. and we will not repeat 
them here. We only note that the continuity equation takes its standard form 

f^+3Hip + P)=0 (30) 

in terms of the Hubble parameter in the Jordan frame, H = dhiR/dt. This is obvious since 
in the Jordan frame, the matter content is not directly coupled to cf). On the other hand: 

p = pJA^ , P = PJA^ . (31) 

These expressions relate the energy density and pressure in the Jordan frame, to those in 
the Einstein frame. 

To proceed, we need to express the pressure and energy density in terms of the tempera- 
ture, through the relations obtained in the previous Section. We can do this in either frame. 
Indeed, the pressure and energy densities, as well as the free energy, can be written as 
times a function of T/m (for example, in the SM , Eq. fll8l) . or msusY in the MSSM, 
Eq. f ll9p ). From Eq. (131 p we deduce that the temperatures T in the Jordan frame and 
in the Einstein frame are related to each either by T = T*/y4. An analogous relation takes 
place between the mass of a particle in the two frames, m = mtf/A. ^ Therefore, the ratio 
between the temperature and the mass scale is the same in both frames, and we can use 
relations like f|T8l) and f|T9l) in either frame, with the rescaling of T and rrit (or msusv)- 

For definiteness, we use these relations in the Jordan frame, where the mass scales are 
constant. Inserting the expressions ([7]) for the energy and pressure in the continuity equation 
(EOI), we get 

A dT ( dgt \ d ( dgr \ 3 dR / dgf \ 

Since gf is a. function of the temperature, 

dgf ^ dgf 1_ dT 
dt dlnT T dt ' 
Inserting this into Eq. ( l32l) . we obtain 

dlnT _ ( 12gf + 3-^ \ dlnR _ 1 dlnR 



dlnT 



(33) 



(34) 



9f. 



^ To see this, consider the action of a massive field in the Jordan frame (where it is of the standard form, 
and the mass parameter is m), transform to the Einstein frame, rescale the field so that it is canonically 
normalized in this frame, and observe that the mass parameter of the newly canonically normalized field 
has indeed become Am = . 
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Using equations (fT8l) and (fT9l) we have 



I=l + '-^al + 0{ar) , SM 

J = 1 + ^^2^0(55/2) , MSSM. 
61 ^ ^ 

In both cases, T — 1 ~ 0.002 at the electroweak scale, and it decreases further as the 
temperature grows (for instance, I —1 = O {10~^) at T = lO^^GeV). If we ignore this 
small effect, and simply set X = 1 , Eq. (!34|) gives us the standard relation T oc 1/R as long 
as there are no other processes involved, such particle annihilations which can reheat the 
thermal bath. Since = AT, and = R/A, the same inverse proportionality oc l/i?* 
also holds in the Einstein frame. 

We can also replace the time derivative on the left hand side of Eq. fl34l) with a derivative 
with respect to p through Eq. (!28|) . Recalling that d\nR/dt = H, and the relation (129|) 
between the Hubble parameters in the two frames, we find 

(,„T)' = -i^ (35) 

where we recall that prime denotes differentiation with respect to p . 

Since a = dlnA/d(f)^:, if we again neglect the small departure of X from one, this relation 
can be integrated to give 

T /2in A (0*,in) _p r T — 1 ^QC^ 

where we recall that, by definition, Pm = . We also have 



, for X= 1. (37) 



The cosmology of the system is therefore completely specified once the evolution (p^ (p) is 
known. Starting from the third expression of (125|) . and replacing the time derivatives with 
derivatives with respect to p through (l28l) . we arrive at 

0: = -^-^ [(1 - 0: + (1 - 3^) a] (38) 



where the equations of state w for the MS and for the MSSM is given in Eqs. ( ITSl) and 
( fT9l) . respectively.^ The equation of state is completely determined by the temperature 
T. Therefore, we can either solve the two equations fl5^ . and fl55]) numerically, or we can 
neglect the small departure of X from unity, and use the analytic expression (136|) for the 
temperature. Our results below are based on the numerical solution, however, we have 
verified that T R = T^R^ remains constant with an extremely good accuracy. 



These are good approximations as long as we do not cross any mass threshold. 
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IV. COSMOLOGICAL SOLUTIONS 



There are strong limits today on deviations from general relativity [30| ■ The bounds are 
most conveniently imposed on the so called post-Newtonian parameters (see Footnote [2]) , 
which in the context of scalar-tensor theories of gravity, can in turn be related to the present 
values ao and Po of the function ([2]), and its derivative 

^(W-^ (39) 

Limits from the Very Long Baseline Interferometer (3l| and the Cassini spacecraft [s^] 
enforce ^ 10~^ . The perihelion shift of Mercury [33| and the Lunar Laser Ranging 



experiment [S^l instead constrain the combination |(/?o + l)«ol ~ 0(10"^). Therefore, j3o 
can be quite large, provided that ao is sufficiently small. As argued in Ref. one should 
however assume that /?o ~ (100), so that the post-Newtonian approximation scheme makes 
sense. 

These stringent limits do not need to hold at earlier cosmological eras, and, in fact, for 
very large /?, a could have been large in the past as well. This is due to the attraction 
towards GR that this type of scalar-tensor theory of gravity possesses This can be 
easily seen from the evolution equation (155]1 for the scalar field. This is an equation for a 
relativistic particle, subject to friction and evolving in a potential 

KfT = ^ (1 - 3«7) 1 a = ^ (1 - 3w) In A (0,) (40) 

(up to an irrelevant constant). Following the conventions in the literature, we define 

A (0.) ^ e'^^'^*) ^ a = -^ , /3 = ^ , Kff = ^(l-3^)a(0.) . (41) 

The evolution naturally tends toward a minimum of the effective potential, for which a = , 
if such a minimum exists. This is more effective for large curvature, namely for large values 
of /9 . Clearly, we need to assume that a has a minimum; this excludes the simplest scalar- 
tensor theory, namely the Brans-Dicke one, for which a is linear in cj)^. In our concrete 
examples below, we consider the simplest example which satisfies the attraction towards 
GR mechanism, namely the quadratic coupling function^ 

a(0*) = f0^ (42) 

The force that drives a towards zero requires that matter fields are coupled to 0* through a 
nonvanishing trace of the stress energy tensor (notice that, indeed, Ke vanishes for w = 1/3). 
The present literature disregards the effect of the trace anomaly pointed out in Section [Tll 
but studies the departure from w = l/?> taking place during the matter dominated era, and 
when the temperature crosses any mass threshold. 



^ In practice, we are Taylor expanding the function a {(j)*) around its minimum, assuming that higher order 
terms are irrelevant. A constant term in a is irrelevant, since it simply changes the constant G* and the 
normalization of the scalar field [cf. the action ((20)) ]. A linear term can be reabsorbed with a shift of . 



12 



It is easy to use Eq. (138|) to estimate the decrease of a during matter domination 
Assuming that the scalar field is sufficiently small so that 0'^ -C a/3 , ^ and that the difference 
between the Jordan and Einstein frames can be ignored in this estimate, Eq. (!38|) is solved 
by a simple damped oscillatory solution. Looking only at the decrease of the amplitude of 
the oscillations, we find 
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g-A (PO-Pcq) 






'^eq 





"eq 



A = Re<^ — 



1/2 ~ 



(43) 



(the suffixes "eq" and "0" denote the values of a quantity at the matter-radiation equality 
era and today, respectively). For > 3/8 (we assume this to be the case in this work), we 
find that a decreases by a factor of about 5 x 10^ during the matter dominated era. This 
estimate can be compared to Figure 11 in the first reference in Ref. l3 ■ 

The decrease when the temperature drops below some mass threshold is more compli- 
cated. Refs. j3| and [l2] studied this effect, in relation to bounds imposed on these theory 
from BBN. Fig. 5 of Ref. [13] shows the "source term" (the coefficient that multiplies 
in the effective potential) due to the mass threshold for any SM particle. One finds one dis- 
tinct peak corresponding to the electron mass threshold, and a series of overlapping peaks 
corresponding to the other mass thresholds, so that there are two distinct phases in which 
0* decreases. The actual decrease of a is strongly sensitive to the values of /?, and 0*, and 



we refer to [l2| for more details. It is possible that a decreases of a factor of O (10~^ — 10~^) 
in either phase. 

We now turn to the departure from w = 1/3 pointed out in Section [TTl We focus on the 
early cosmology before BBN, at temperatures higher than any mass threshold. Specifically, 
we consider temperatures greater than 200 GeV for the SM, and greater than 500 GeV for 
the MSSM (as we mentioned in Section [III we assume that this is the mass of all of the 
supersymmetric partners of the SM).'' 

We show in Figure [2] the decrease of a (0*) due to the trace anomaly. For illustrative 
purposes, we choose the initial temperature Tin = 10^ GeV in the Jordan frame. We then 
evolve the two equations fl35l) and fl38|) starting from several different values of 0* and /? 
using 0'^ in = ~ (9/3£/2) 0*in , according to the slow roll solution of its equation of motion (see 
the discussion in the paragraph above Eq. flHT]) ). The initial value of 0*, and determine 
the initial value ain. This, in turns, gives us the initial value of the temperature in the 
Einstein frame, T*in = exp (oin) Tin. We chose values of a^n ^ 32 , to guarantee that the 
initial temperature in the Einstein frame is sub-Planckian. In each case, the evolution is 
concluded when the temperature in the Jordan frame reaches 200 GeV (500 GeV) in the SM 
(MSSM) case, so that we do cross any mass threshold. We denote by Oout the value of a at 
this moment, and we plot in Figure [2] the ratio ctout/am vs the initial value Oin. 

The ratio ctout/flin is nearly constant for small values of Oin. The reason is the following: 
small values of a^^ correspond to small 0,,. In this case, Eq. ( |38|) for 0^, is approximately 



^ In this case, Eq. (j38|) is a linear equation, and the ratio between the final and initial values of 0* is 

independent of the initial value of 0* . 
^ The numerical values refer to the temperature in the Jordan frame, since this is the frame in which the 

particles have constant mass; notice however that the ratio between the temperature and any mass scale 

is identical in the two frames, see the discussion after Eq. (I^T]) . 
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FIG. 2: Left panel: The decrease of a = (3<p1/2 for the SM, starting from different values of 
a at the initial temperature T = 10^ GeV in the Jordan frame, and ending at the temperature 
T = 200 GeV . Right panel: Decrease of a in the MSSM with the same initial temperature, but 
with the final temperature T = 500 GeV. See the main text for details. 



linear, since one can neglect 0'^ at the denominator 
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(44) 



so that the ratio out/0* in is independent of the initial normalization of 0* . Since a oc (pl, 
this can be rephrased in the statement that aout/om is independent of the initial value of 
ttin . We can also use the approximate equation (jHj) to understand why the decrease of 0* 
observed in Figure [2] is so sensitive to the value of /? . 

From Equation (jS]), we see that if e vanishes, the scalar field freezes to a constant value 
typically within one e-fold of expansion. We are thus in a over-damped regime in which the 
velocity is almost constant and fixed by the source term (see below for further justification). 
Then, as long as the field is slow-rolling and the deviation of T from 1 is small, Eq. 
leads to 

9 



-pedlnT. 



(45) 



Neglecting the variation of e in the l.h.s. of fHSl) (since it is suppressed with respect to the 
variation in the r.h.s.), one can easily integrate this equation to get 
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Equation (H6l) can be solved in terms of the Lambert W function, that is the inverse function 
of w — s> wexpw, as 
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(48) 
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where we have chosen to fix e to its geometric mean e = y e{Ti^it)e(Tout) in the l.h.s. of 
Eq. ( H5l) before integrating. This expression reproduces the computation depicted in Fig. [2] 
with great accuracy, for P < 50, and < 10 . 

Another way to grasp some insight on the evolution of the scalar field is to parametrize 
the solutions of Eq. (jS]) as 0* = const, x exp A {p') dp' . The function A must then satisfy 

X' + X'+^l + ^e^X + ^ep-O. (49) 

Let us emphasize that e > in both the SM f |T8l) and the MSSM ( |T9|) . The equation of 
state w is slowly evolving and close to that of radiation (see Figure [1]). Therefore, e' and e"^ 
are both much smaller than e. We then find the two solutions: Ai ^ —1 and A2 — —9 (3 e/2 . 
This gives us two approximate solutions of the linearized equation fH4|) . To understand these 
solutions, consider the case e = in which the effective potential for 0^, is absent. The two 
solutions then give 

0, = Ci e/' '^P' + C2 e^" '^P' = Ci e-(f-P'") + C2 = Ci?^ + C2 , £ = 0. (50) 

Namely, starting from generic initial conditions (arbitrary integration constants), quickly 
dissipates its velocity, due to the friction term in its equation of motion, and stops at the 
constant value C2 ■ For small £ , as we are considering here, the change in the fast decaying 
solution can be neglected; however, the constant mode oc C2 is now replaced by a slowly 
decreasing mode. This solution is a slow roll solution of the linearized equation (jS]), in 
the sense that the (f)'l term can be neglected (as well as the term oc £ ^ 1 multiplying cf)'^ ). 
It is clear that, starting from a generic initial 0^,, after a quick initial transient regime the 
fast decreasing mode will become negligible. Therefore, in our numerical computations, we 
assume that 0* is in the slowly rolling regime initially. Therefore, as long as Eq. (jUj) is a 
good approximation, 9/?e/2 ^ 1 , and e evolves sufficiently slowly, we have the approximate 
solution 

0,(p)^0,,i„e-f^r*'^(p'). (51) 

We see that the final value of 0* is indeed exponentially sensitive to (3 . 

The accuracy of the approximate solution fl5Tl) is shown in Figure [3] for two specific choices 
of the parameters. We take sufficiently small so that the linearized Eq. flH|) holds, and we 
choose two different values of (3 (we also choose Tin = 10^ GeV, and we conclude the evolution 
when T = 500 GeV, as in Figure [2]). Accuracy of the approximate solution requires that 
9/2/5 — 0.002 /? <^ 1 . For /3 = 50 , the two curves differ of about 2% at the final moment; 
for P = 100, the final discrepancy is about 20% . 

We also see in Figure [2] that a decreases more if we start at greater a[^. Moreover, the 
decrease is stronger for the MSSM than for the SM. It is easier to understand this effect 
by discussing the evolution in the Jordan frame. As the scalar field evolves towards the 
origin, a(0*), and A = exp a, decrease. If A decreases sufficiently fast, the temperature 
in the Jordan frame can actually increase. To see this, recall the relation between the 
temperatures in the two frames, T = T^,/A , and the fact that the cosmological evolution is 
standard in the Einstein frame. As we discussed at the end of Section llllt the temperature 
in each frame is (up to O (10~^) corrections) inversely proportional to the scale factor in that 
frame. Therefore, if A (0^,) decreases sufficiently fast, the universe contracts in the Jordan 
frame. As the temperature rises and then decreases again, the scalar field experiences a 
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FIG. 3: A comparison of the exact and approximate solution ()5ip for <j)^,. The two sets of curves 
correspond to the same initial temperature and a = (3 4>l,/'2 , but to two different values of (3. 
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FIG. 4: Left panel: The evolution of the temperature in the Jordan and Einstein frame, for a 
specific choice of parameters. The temperature in the Jordan frame is increasing at the earliest 
times shown, it reaches a maximal value, and it then decreases. The universe contracts while 
T increases. Notice also that the two frames coincide at late times. Right panel: The maximal 
temperature reached in the Jordan frame for different values of the parameters (chosen to be the 
same as in Figure [2]) . 



"driving force" towards the origin for a longer period of time. This causes a larger overall 
decrease in (j)^ with respect to the cases in which (p^ (and, hence A) is initially small, and the 
evolution is standard in both frames. The suppression is more marked in the MSSM than 
in the SM, since the equation of state has a greater departure from 1/3 in that case (this 
is particularly true at higher values of T, so at larger A-^^ , cf. Figure [1]) . ® A comparison 



As we already mentioned in Section IlIIl the key quantity determining the departure from w = 1/3 is 
actually the ratio between the temperature and the particle mass scale, which is the same in both frames; 
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between the right panels of Figures [2] and H] shows that indeed the ratio ctout/flin starts 
decreasing, when ain starts out sufficiently large so that the Jordan frame temperature is 
initially increasing. 

In the left panel of Figure H] we show the evolution of the temperature in the two frames 
for some specific choice of the parameters. We choose and (3 sufficiently large, so that the 
temperature in the Jordan frame increases for some time. We verified that the temperatures 
are inversely proportional to their respective scale factors in each frame with a very good 
accuracy (more quantitatively, we found that the product T R = T^, increases by about 
1% between the initial and final times shown). This is the reason why log^g vs. p oc log i?* 
appears as a straight line in the Figure. We know from Figure [2] that a experiences a strong 
decrease for this choice of parameters. Indeed, a ~ 0.003 at the end of the evolution; 
therefore, the two frames, and the respective temperature, are (nearly) coincident at late 
times. 

In the right panel of Figure H] we show instead the maximal temperature T reached 
for different choices of parameters. In general, the temperature reaches a higher value for 
greater values of and /5 . These are the cases for which a shows a strong decrease, cf. the 
right panel of Figure [2l It is easy to estimate when the temperature in the Jordan frame 
increases. From Eq. (!35|) . and recalling that X is always very close to one, we have that the 
temperature increases whenever 0*0*' < —l/P- Using the slow roll relation 0'^ ~ —9eP/2 , 
and our choice (142|) of the function a, this condition becomes 

a > — — =^ T increases . (52) 
Q p £ 

We can compare this estimate with the result shown in the right panel of Figure |H Using 
T = 10^ GeV in Eq. f lT9|) . we have an initial value of e-m ^ 4.2 x 10""^ , and we find that 
the temperature is initially increasing provided that 270/(3. More specifically, the 

Jordan frame temperature is initially increasing when ^ 270, 27, 5.4, 2.7, 1.8 , for (3 = 
1, 10, 50, 100, 150 , respectively, and hence T^ax will be greater than the initial temperature 
of T = 10^ GeV. This is very well con&med in the right panel of Figure |H 

It is clear that the Jordan frame temperature always reaches a maximal value, and then 
decreases. Indeed, even assuming that the inequality flS^ is initially satisfied, a oc 0^ 
decreases, while the right hand side slowly increases (since e slowly decreases as T grows). 
Tmax is reached whenever the inequahty stops holding. Therefore, a\T^^^ — 1/ (9/?e), and 

0*1^^^^ ~ \j2/e/ (3 (3) when the Jordan frame temperature reaches its maximal value. If we 
also have an analytical expression for p at this moment, we can then use Eq. (!36l) to obtain 
an analytical expression for T^ax • The redshift p can be related to the value of 0* that 
we have just obtained through the approximate solution (15T]) . Unfortunately, e varies with 
temperature, so that we cannot straightforwardly invert Eq. (!5T|) to find p as a function of 
0*. To obtain an estimate, we simply assume that e is constant, so that 

0* \^ , , 
, e = const. (53) 

0*,in / 

By inserting our estimate for 0*|Tmax) ^ind the relation between 0*^in and in this expression. 



so, this discussion can be equivalently done in the Einstein frame. 
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we obtain an analytic expression for p at Tmax- Inserting all of this into Eq. (l36i) . we find 



This estimate assumes that the inequality fl52l) is initially satisfied, so that T is growing 
initially (if this is not the case, Tmax is simply Tin). A comparison with the exact values 
for Tmax shown in the right panel of Figure H] indicates that the estimate fl5^ is accurate 
provided T^ax is about 2 — 4 orders of magnitude greater than Tin (the precise value depends 
on the choice of (3). For a larger span of temperatures, e can no longer be assumed to be 
constant, and the estimate fl5^ loses its accuracy. 

We conclude this Section with a few comments on the effect of starting from an arbitrary 
value of (f)'^ (rather than starting from the slow roll value obtained from Eq. flHTl) ). One way 
to ensure that 0* is always in the slow roll regime is to imagine that the field responsible 
for inflation is coupled differently to 0* than to matter fields. For instance, if the infiaton 
is decoupled from (p^, in the Einstein frame, 0^, has no driving force during inflation. Due to 
friction, 0^. reaches a constant value, and inflation proceeds as in GR. The effective potential 
for 0* is generated at reheating, when the infiaton decays to the matter fields which are 
coupled to the scalar. In this scenario, it is reasonable to assume that 0* will settle in the 
slow roll regime as reheating proceeds and the effective potential is being created. However, 
without a complete theory, it is legitimate to assume any value for 0'^ at the start of our 
simulations. 

Wen considering arbitrary initial velocity, we have however to restrict the choice to |0'^| < 
V^. Indeed, the first equation of fl25l) can be rewritten as: i?^(3 — (p'^) = SttG^p^. Solutions 
with 10'^ I > require p^, < and are therefore unphysical. Moreover, the two branches 
of solutions with |0'^| greater or smaller than \/3 are mathematically disconnected. For the 
physically meaningful case, the effect of the 3 — (j)'^ prefactor in Eq. fl38l) can be neglected 
in the present qualitative discussion, so that with our choice of a, we simply have the 
equation of a scalar in a quadratic potential subject to friction. The friction will damp the 
initial velocity, so that, after an initial transient regime, the field proceeds as in the cases 
we studied. We can however imagine the situation in which in this initial phase the field 
is actually moving away from the origin (it will then stop, due to the friction term, and 
move back towards the minimum at 0* = 0). In general, as is it clear from Eq. (136|) . the 
temperature T in the Jordan frame decreases faster as |0=k| increases, while it decreases more 
slowly (or even increases) when |0*| decreases. 

V. DISCUSSION 

In this work, we have studied the dynamics of CCSFs in a radiation dominated universe. 
We have shown that there are cases in which the evolution of the scalar field can have a 
very strong impact on cosmology. Despite the conformal coupling, the scalar field has a 
source term given by the conformal anomaly, which is ultimately related to the fact that 
particles in the thermal bath are interacting (this leads to running coupling constants, and 
to T^ = p — 3 P 7^ in the thermal bath) . We computed the leading value for the source 
term both in the SM and the MSSM. We then focused on the effect that it can have in the 
context of scalar tensor theories of gravity, in which the scalar field (p^ is conformally coupled 
to matter (by "matter", we mean any SM or MSSM field). We specifically considered scalar 




e = const. 



(54) 
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tensor theories which have GR as an attractor. In such theories, the additional couphng 
between matter mediated by 0* has a strength which is function of the scalar field itself, 
and which vanishes for some given value of 0*. One simple way to formulate the theory 
is to impose a standard (^^-independent) action for the matter fields, but a nonstandard 
curvature term, with oc (0) - this is known as the Jordan frame. Then, the strength 
of the new interaction is controlled by a = d\nA/d(f)^. If present, matter fields with 7^ 
tend to drive (p^, to a value for which a vanishes, and GR is recovered The efficiency of 
this mechanism increases with increasing T^, and with increasing curvature of In A close to 
the minimum. The latter quantity is typically denoted by (3 ; our explicit study is performed 
for the simplest case in which j3 is constant (^^-independent). However, we expect that the 
qualitative features that we find are general. As we remarked, precision tests of GR impose 
that al < 10~^ and \Po a^l < O (10~^) today. This does not implies that P should be small 
(indeed, increasing f3, leads to a smaller value for a today). Moreover, this leaves significant 
room for departures from GR at earlier cosmological eras. 

As we showed in Figure [T], the trace anomaly gives a small departure from the noninter- 
acting result, 1/3 — w = O (10~^ — 10~^), at temperatures much greater than the mass of 
any SM or MSSM particle. Nonetheless, in some cases this can have a strong impact on the 
evolution of 0*, and on the whole cosmology. The most obvious effect is that a decreases 
during this era, due to the mechanism of attraction towards GR [3]. The amount of the 
decrease is depicted in Figure [21 where we see that the effect is strongly sensitive to j3 and 
to the initial value of A. Previous cosmological studies of such theories (see for instance 
Ref. [nl, T2|) only considered the attraction towards GR during (i) the matter dominated 
era {w = 0), or (ii) in the radiation dominated era, whenever the temperature crosses any 
mass threshold (this results in w 7^ 1/3 when the temperature is close to the mass of any 
particle in the thermal bath). Our work refers to the previous cosmological epoch, and 
therefore, our final values should be understood as the initial values for such studies. Since 
a in general decreases for 7^ 0, we generally extend the region of allowed initial conditions 
for such models. The amount of the extension depends on the initial temperature, as well 
as on the initial values of (3 and a . To provide some "benchmark" , based on the example 
we studied, we found that our effect is non-negligible whenever f3 10, and very strong for 
/3 ~ 50 — 100 or greater. 

A second, less obvious, effect of the conformal anomaly is that it can lead to a contract- 
ing phase of the universe, followed by a regular bounce, and standard expansion. Similar 
(bouncing) behaviour has been postulated as a possible alternative to inflation and in other 
cosmological contexts as well [37|. The cosmological evolution appears as standard in the 
Einstein frame (the frame in which the gravity term is regular, but the metric felt by mat- 
ter contains the function A), with 0* contributing to the energy density and pressure as a 
standard scalar fleld. Therefore, the scale factor of the universe R^, in the Einstein frame 
always increases. The temperature T in the Jordan frame is related to the temperature T* 
of the Einstein frame by T = T^,/y4 (0*). Moreover, from the conformal transformation we 
have R = R^: A (0*) . If the evolution of 0^, leads to a sufficiently fast decrease of A, then 
R can decrease, even if R^: is increasing. This leads to a temporary phase in which the 
universe contracts, and T increases. Also in this case, the effect is more marked for large (5 
and = e""\ This can be seen in Figure HI where we show the maximal value reached by 
T when the universe stops contracting and bounces, for different values of (3 and Ain, but 
for the same initial temperature T = 10^ GeV. We see that T can increase by several orders 
of magnitude before the bounce. 
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Before discussing the possible consequence of the bounce, it is useful to make a few 
clarifications. The first is that whether the universe is contracting or expanding is obviously a 
frame independent question, whose resolution is most immediately seen from the behavior of 
R rather than R^,. To see whether the universe expands or contracts, cosmological distances 
need to be compared with a "ruler". It is convenient to choose the Compton wavelength 
of any particle, A oc m"^ (where m is the mass) as the ruler j36|. In the Jordan frame, 
A is constant (since the action of matter is independent of 0^,). Therefore, the universe 
expands/contracts whenever R increases/decreases. It is easy to show that the mass of a 
particle m* in the Einstein frame is related to the mass m in the Jordan frame by m* = 
mA (0*) . Therefore, the ratio between the scale factor and the Compton wavelength is the 
same in both frames. The second clarification is that, in principle, one can have a regular 
bounce in a scalar tensor theory of gravity even if one disregards the trace anomaly. Since 
T = T^:/A (0*), one "simply" needs to arrange for a sufficiently fast decrease of A for some 
period of time. For instance, it may be possible to obtain this with a sufficiently high initial 
velocity for (j)^ . Here, we have neglected any transient phase due to the initial value of 0'^ . 
From the evolution equation for 0* one can easily see that the initial velocity is dissipated 
in about one e-fold of expansion (in the Einstein frame). After that, the evolution of 0* is 
simply dictated by its coupling to matter. In the early universe, neglecting the conformal 
anomaly leads to 0'^ = . Including (as one must) the effect of the conformal anomaly leads 
to the bounce studied here. 

Regarding the possible consequences of the bounce, it is tempting to wonder whether it 
could lead to a solution of some of the standard cosmological problems, such as the horizon 
problem. If the universe could contract and remain small for a long enough time, one could 
hope that the propagation of light in that time could result in a much larger horizon than 
achieved in a standard cosmology. We have verified that this is not the case. This can 
be most easily understood in the Einstein frame, in which the cosmological evolution is 
standard, leading to a horizon size which is nearly coincident with the standard result. We 
can rephrase this by stating that, in the Jordan frame, the universe remains small for a 
negligible amount of time. 

The fact that the universe bounces and reaches a maximal temperature in the Jordan 
frame can have instead very interesting consequences for the thermal production of parti- 
cles and topological defects. Particle physics processes are most immediately computed in 
the Jordan frame, in which the particles masses and interactions are independent of 0^, . 
Therefore, computations of thermal scatterings and symmetry restorations are standard in 
the Jordan frame. We can arrange T^ax to remain sufficiently small, so that the thermal 
production of unwanted relics is sufficiently suppressed, and no symmetries which would 
result in unwanted defects are restored. This could happen even if the temperature 
in the Einstein frame is high, see for instance the left panel of Figure HI Provided is 
sufficiently high, so to provide an initial large hierarchy between the Einstein and Jordan 
frame temperatures, and (3 is not too large, so that T does not increase too much during the 
evolution (see the right panel of Figure H]), one could even start with the highest possible 
value ~ G* in field theory (i.e., before entering in the quantum gravity regime), with- 
out T ever conflicting with any cosmological bound, such as that arising from the gravitino 
problem. The presence of a bounce could also have consequences for wanted relics, namely 
baryons and dark matter, if Tmax is close to the decoupling temperature of dark matter or 
to the mass of a particle responsible for the baryon asymmetry (for instance, the lightest 
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right handed neutrino in leptogenesis) . ^ 

We stress that, while the present study has been focused on scalar tensor theories of 
gravity, the effect we have studied is relevant for any CCSF. Another interesting application 
is the evolution of the dilaton $ in string models. The dilaton is also conformally coupled 
to matter, and it posses a mechanism analogous to the attraction to GR we have studied, 
namely the least coupling principle [sj. Also in that context, one can go to the Einstein 
frame, in which the gravity term is standard, and canonically normalize the matter fields. 
Neglecting the interactions between different fields, the dilaton then enters in the expression 
for the masses of the various fields. This results in the following equation of motion for the 
dilaton ^ 



3 - 

$" = _r — 

2 



[1-w) ^' + J2aA 



Pa-SPa 



«A = (55) 

where prime again denotes derivative with respect to p, the index A denotes any matter 
field, and p is the total matter energy density, p = J2aPa- This equation is analogous to eq. 
f l38|) for 0*. If there was only one matter field, the expression would actually be identical. 
Indeed, since in scalar tensor theories of gravity = A{(j)^,) m, with a 0^,— independent 
Jordan frame mass m, the coupling a entering in Eq. (155]) and defined in Eq. ([2]) can also 
be expressed as a = d In m^/dcl):^ . For the dilaton case, there is the additional comphcation 
that the functional dependence {^) can be different for different matter fields. However, 
if all the functions aA vanish for one specific value of then the matter fields drive the 
dilaton to that value, in the same way as 0* is driven to a zero of a jsj. For this value, the 
dilaton decouples from matter. Due to the trace anomaly we have discussed, the dilaton 
will be attracted to that value also during the radiation dominated era. This can possibly 
result in larger range of initially allowed values for $ than what one would have neglecting 
the trace anomaly (as for the case we studied, the size of this effect depends on the specific 
functional forms of ($) ). 

As we mentioned, the consequences of the conformal anomaly are more sizeable for large (3 
and y4in. In most of the examples shown, we started from Tin = 10^ GeV, and we considered 
y4in up to ~ 10^^ , so that, when this limit is saturated, T*^in is Planckian. It is tempting to 
relate the hierarchy between the electroweak and the Planckian scale to a large conformal 
factor. This is precisely what is done for the dilaton in the context of the so called little 
string theory at the TeV [s^. Consider the string theory expression for the Planck mass 
Mp ~ M^Ve/g^, where gg is the string coupling, the string scale, and Vq the volume 
of the six-dimensional internal space. Assuming that the fundamental scale is TeV, a large 
4d planck mass can be obtained by a large volume Vq 13]. Alternatively, one can assume 



that both Ms and Vq are set by the TeV scale, and that the large Planck mass is due to a 
very small string coupling gg ~ 10^^^ . From the expansion series gg = e'^ + . . ., this then 
requires a large negative value for the dilaton. Therefore, it is argued in the last of Ref. 
that the large value of the Planck mass today can be related to the runaway of the dilaton, 
assuming that nonperturbative terms in the potential stabilize $ at the desired value (this 



9 



Limits on scalar tensor theories from the dark matter relic abundance were obtained for instance in 
Ref. 



38[ by considering only the evolution of 0, due to the mass thresholds. 
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work also discusses limits on the theories imposed by deviations from GR). In our context, 
we can express the Planck mass (i.e., the coefficient of R in the action) in the Jordan/string 
frame as (cf. Eq. (|20|)) 



2 



= - - = — = TeV^ K^^-'f'') (56) 

Namely, we assume that the Planck mass is initially at the TeV scale, and it is then driven 
to a large value by the evolution of 0*. Therefore, the last term in fl56l) plays the same 
role as I/qs ($) in the string theory expression. Although our choice of A (0*) was dictated 
by simplicity, one interesting extension of our study would be to consider other functional 
forms which reproduce the string dilaton case. 

To conclude, modifications of GR can result in interesting nonstandard cosmologies at 
early times, when such modifications are less constrained. In this work, we considered the 
effect of the trace anomaly at temperatures higher than the electroweak scale. Somewhat 
surprisingly, we saw that that the universe can experience a regular bounce in this context. 
As we approach this early time regime, it is natural to study how this effect can be reconciled 
within a complete cosmological framework, where the standard cosmological problems are 
solved, and the primordial perturbations are obtained in agreement with observation. We 
believe that this issue deserves further investigation. 
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